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^ ■ Abstract 

O. 

We present the complete calculations of the forward-backward asymmetry (^fb) and the total 



(N 



cross section of top quark pair production induced by dimension-six four quark operators at the 



i-C ' Tevatron up to 0(a^/A^). Our results show that next-to- leading order (NLO) QCD corrections 

o ■ can change ^PB and the total cross section by about 10%. Moreover, NLO QCD corrections 

,^ , reduce the dependence of ^fb and total cross section on the renormalization and factorization 

_i ^ scales significantly. We also evaluate the total cross section and the charge asymmetry (Aq) 

pi ■ induced by these operators at the Large Hadron Collider (LHC) up to ©(a^/A^), for the parameter 

T-H ! 

(^ . space allowed by the Tevatron data. We find that the value of Ac induced by these operators is 

I I , much larger than SM prediction, and LHC has potential to discover these NP effects when the 

measurement precision increases. 
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I. INTRODUCTION 



The top quark is the heaviest particle discovered so far, with a mass close to the elec- 
troweak symmetry breaking scale. Thus it is a wonderful probe for the electroweak breaking 
mechanism and new physics (NP) beyond the standard model (SM) through its productions 
and decays at colliders. The forward-backward asymmetry (Afb) of the top quark pair 
production is one of the interesting observables at hadron colliders. Within the SM, Ap-Q 
is absent at the tree level in QCD due to charge symmetry, and occurs at next-to-leading 



M- In 
lo|. Re- 



order (NLO) 0{al) in QCD with the prediction ApB ~ 6% in the ti rest frame 
the last few years, D0 and CDF Collaborations measured Ap-Q at the Tevatron |7|- 
cently, the CDF Collaborations annouced that, for the invariant mass of the top quark pair 
TTT'tt > 450 GeV, the measured asymmetry, Afb = 0.475 ± 0.114J9|, differs by 3Aa from the 
SM predictions A-pB = 0.088 ± 0.013, which has aroused many discussions of explaining this 
deviation in NP model, including new gauge bosons, axigluons and so on ll|-l68|. 

Since we do not know which type of NP will be responsible for this deviation, it is interest- 
ing to study the A-p-Q in a model independent way, using an effective Lagrangian. In general, 
NP scale relevant to Ape is large enough so that the heavy fields have been integrated out 
at the low energy scale. At the Tevatron, the subprocess qq — )■ ti dominates over top quark 
pair production, so only contributions from dimension-six four quark operators to the ti 
production are considered. Similar approach had been adopted for the dijet production to 



constrain the composite scale of light quarks [69 

be written as 



731]. The relevant effective Lagrangian can 



^^P = J^Y1 [^AB {ui.qA) iiBl^ts) + C% iuT'^l^qA) iisT'^rtB 



:i) 



A,B 



where {A,B} = {L,R} with q = {u,d)'^, (c, s)^. Up to C(as/A^), the NP contributions to 



the total cross section and the Ap-Q are clear in the vector-axial basis 28 



,Q, 



as compared 



with the chirality basis. Only the axial- axial current combination contributions to the Ap^, 
and the vector-vector operator contributes to the total cross section. However, this is no 
longer true up to 0{a^/A^). The chirality basis is the preferred one when studying the 
chiral structure of NP effects much above the Electroweak scale, so we choose to work in the 
chirality basis. The contributions to Arb at leading order (LO) from such operators has been 



explored in Refs. 



19 



28 



29 



31 



45l . |75| . It is shown that the Ap^ observed at the Tevatron 



can be explained by above operators for suitable parameters. As we know, the LO cross 
section at hardron colliders suffers from large uncertainties due to the arbitrary choice of the 
renormalization scale and factorization scale, thus it is important to include NLO corrections 
to improve theoretical predictions. Besides, at the NLO level, virtual corrections, real gluon 



ead to a sizeable difference between the 
2|, which will contribute to ApB- 



emission and massless (anti) quark emission can 
differential top and anti-top production process I 

In this paper, we present the complete NLO QCD calculations of A-pB and the total cross 
section of top quark pair production at the Tevatron induced by above operators, and we 
also study the top quark pair production at the Large Hadron Collider (LHC) induced by 
these operators at the NLO QCD level. Last year, LHC reported their first observation of 
top quark pair production, and will soon become a major top quark factory. At the LHC, the 
top quark pairs can be produced through quark antiquark annihilation qq — )■ ti and gluon 
fusion gg — )■ ti. Since gluon fusion channel dominates at the LHC, it is difficult to probe 
these four quark effective operators from early LHC results. However, it is still possible 
to detect these effects from above effective operators on the Charge Asymmetry (Ac) at the 
LHC, in the parameter space allowed by the Tevatron data, when the measurement precision 
increases. 

The arrangement of this paper is as follows. In Sec. Ullwe show the LO results. In Sec. IIIIl 
we present the details of the NLO calculations, including the virtual and real corrections to 
the top quark pair production. Section IIVI contains the numerical results, and Section |V] is 
a brief summary. 

II. LO RESULTS 



Throughout our calculation, we adopt the same conventions as in Ref. 76| (see Sec. HI 
A), and present the he 
regularization scheme 



icity amplitudes for qq — )■ tt in the Four-Dimensional Helicity (FDH) 
77| ■ The ti production amplitudes, including NP contributions, can 



be written as 



Mtt = asko + ^/lo + "s/nlo + ]y^/nlo + " " " ' (2) 



and thus the partonic cross section, up to 0{a1/ h?), can be written as 



a. 



- ^2/SM/.SM* I o"'5'0^ ^-fSM A-NP*\ 

^tt = "s/lo/lo + 2^/^e (/lo /lo j 



+2alne [fl^f^Z) + 2g [JZe {f^^ f^^o) + ^e {fl^f^lo)] ■ (3) 

The LO Feynman diagrams for the subprocess q{pi)q{p2) — ^ '^(^3)^(^4) induced by the SM 
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FIG. 1: LO Feynman diagrams for qq — t- tt induced by SM QCD and NP interactions. 



QCD and the NP interactions are shown in Fig. [H and their (H 1-+) hehcity amphtudes 

are 






++) = '-^iMi + M2)Cs, 



s 
i2mt 



(4) 



[(MiC^R, + A^^CAl) Ci + {M^Cl^ + M^Cl^) Cs] , (5) 



where the SM QCD and NP contributions are denoted by superscipts SM and NP, and 
Mandelstam variables s, t and u are defined as follows: 



s = ipi+ P2)^ t = {pi~ P3)^ u = (pi- PiY- 



We define the following abbreviations for the color structures and matrix elements. 



(6) 



Ml 



TWc 



(^4l)(r/3|3|2] 

fa2)(773|3|l] 
Ci = Si^i^d. 



M, 



Ma 



(^3l)(r?4|4|2] 

(3^773) (r/44^)^ 

(^32)(r/4|4|l] 
(3^r73)(r/44^)^ 



«3«4! 



C.=TLTL 



«2«1 «3«4' 



(7) 



where ii,..4 are the color indices of the external quarks and the boldface momenta denotes 
massive vectors. We use the modified spinor helicity method suited for massive particles |78| 
in our calculations, and a recent application of this method can be found in the Ref |79| . 



The ( — h ++) amplitudes are given by 

Ml^i- + ++) = ^^^ (M, + M,) Cs, (8) 

-^Lo(- + ++) = ^ [{MsCU + M,Cl^) C, + {MsCl^ + M,Cl^) Cg] . (9) 

The amphtudes with other hehcity configurations can be obtained from (H (-+) and 



' — h ++) by exchanging hght-hke momenta p and rj [76|, |79[. At the LO, there is only 



m. 



vector current coupling ip^'^tp at the massive quark vertex. At the NLO, however, magnetic- 
momentum coupling ilj{ia^'^{p3 + p^) ^)ip / {2mt) is induced from loop diagrams. For com- 
pleteness we list the matrix elements for magnetic-moment interaction, 

(„) ^ m,^(r/3l)(r74l)[21] („) ^ (12)(r^3|3|2](r^4|4|2] 

(3^r^3)(r/44^) ' ' (3^r/3)(r744^) ' 

Am)_ ml{r^,2){T^,2)[2l] (^) _ (12)(r/3|3|l](r/4|4|l] 

After phase space integration, the 0{as/K^) partonic differential cross section is 



dcosO ~ 18 A2 



1(1 + p + p^cos'e)icl^ + cl^) + l^coseicl^ - cl^) 



:ii) 



where p = Amf/s, (3 = a/1 — p^, and 6 is the polar angle between the incoming quark and the 
outgoing top quark in the ti rest frame. The color and spin indices are averaged(summed) 
over initial(final) states. In Eq. ( ITTj) the term linear in cos 6 could generate Ap^ proportional 
to (C|^j^ — Clj^) and the rest terms contribute to the total cross section proportional to 
(^RR + ^lr)- These relations will be changed at the NLO level. 

The LO total cross section at the hadron collider is obtained by convoluting the partonic 
level cross section with the Parton Distribution Function (PDF) fi/^ for the initial hadron 
A: 

O-LO = X] / dXa dXbfa/A{Xa,fJ'f)fb/B{Xb,fJ'f)crLO, (12) 

a,b -^^ -^^/^-^ 

where r = 4m^/s. The sum is over all possible initial partons. 

III. NLO QCD CORRECTIONS 

The NLO corrections to the top pair production consist of the virtual corrections, gen- 
erated by loop diagrams of colored particles, and the real corrections with the radiation of 



a real gluon or a niassless (anti) quark. We carried out all the calculations in the 't Hooft 
Feynman gauge and used the FDH scheme to regularize all the divergences. Moveover, for 
the real corrections, we used the dipole substraction method with massive partons 



to 



separate the infrared (IR) divergences, which is convenient for the case of massive Feynman 
diagrams and provides better numerical accuracy. 

A. Virtual corrections 



The virtual corrections for the top quark pair production include the box diagrams, 
triangle diagrams, and self-energy diagrams in SM QCD and NP as shown in Fig. [2] and 
Fig. [31 We have calculated the one- loop helicity amplitudes for the SM process, and find 
complete agreement with those in the Ref . 76|, |8l| ■ Here we only list the NP contributions. 
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FIG. 2: One-loop virtual Feynman diagrams for qq — )• tt induced by SM QCD interactions. 



All the ultraviolet (UV) divergences in the loop diagrams are canceled by counterterms 
for the wave functions of the external fields {6Zq,6Zt), and the Wilsion coefficients 6Z(ji . 
For the external fields, we fix all the renormalization constants using on-shell subtraction. 
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FIG. 3: One-loop virtual Feynman diagrams for qq — )■ tt induced by NP interactions. 



and, therefore, they also have IR singularities 
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-31n^ 



(13) 
(14) 



where C^ = {Any ^,^_ . , Cp = 4/3 and fir is the renormalization scale. For contourterms of 
the Wilsion coefficients SZ^n , we adopted the MS scheme 
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(15) 



/ 



where n/ = 5 is the number of massless quarks appearing in the closed loop diagram, and 
the order of the Wilsion coefficients is 



//-il ^8 r^\ ^8 ^1 ^8 ^1 ^8 \ 
y^lA,! '-'LL; '-'LR; '-'LR) '-'RL; '-'RL; '-'RR; '-'RR^ • 



(16) 



We have considered mixing effects of different color and chiral operators, and the evolution 
equations of the Wilson coefficients are given in the Appendix |Bl The renormalized virtual 
amplitudes can be written as 



M^ = M' 



+ M' 



(17) 



Here y\4'^°''^° contains the self-energy and vertex corrections, and A^^°" are the correspond- 
ing counterterms. The renormalized amplitude A^^ is UV finite, but still contains IR 
divergences, which are given by 






(18) 
(19) 



where we define the IR divergence coefficients /gf", f^^, flf and f^ for different color 



configurations, "s" for singlet and "o" for octect, 
1 1 



IR 

ss 






^IR 



iH>m 



+ ln 



/i: 



IR 

so 



IR 

OS 



2 \uij eiR 

3 \uij eiR 



(20) 
(21) 
(22) 
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oo 
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^IR ^IR 
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/3 



^mf^Ulin^"^ 



IV fP + i 



In 






In 



fij 



4 V/^r/ 2 \fi^ 

where ti = mf — t, Ui = mf — u and A^f^, M^^ and M.^ are defined as follows. 



(23) 



<NP 



NP/ 



MZ = M^Ci + M^rCs. 



(24) 
(25) 



Since we only consider high order corrections up to 0{a^/A'^), the IR divergences of the 
virtual corrections can be written as 

27^e [yW^sM-^N?1 + 27^e [M'^pM^^*] 
'-C,C, [{9C + 2flf) 7^e (Air*-Mf ) + 4/^^7^e {MrMf")] . (26) 

The finite terms in A^^p ^^^ given in the Appendix |X1 



TT 



B. Real corrections 



At the NLO level the real corrections consist of the radiations of an additional gluon or 
massless (anti) quark in the final states, including the subprocess 



qq -^ ttg, gq{q) -^ ttq{q) 



(27) 



as shown in FigJUand FigJS) 












FIG. 4: Feynnian diagrams for the real gluon emission contributions induced by SM QCD and NP 
interactions. 



Before performing the numerical calculations, we need to extract the IR divergences in 
the real corrections. In the dipole formalism this is done by subtracting some dipole terms 
from the real corrections to cancel the singularities and large logarithms exactly, and then 
the real corrections become integrable in four dimensions. On the other hand, these dipole 
subtraction terms are analytically integrable in n dimensions over one-parton subspaces, 
which give e poles that represent the soft and collinear divergences. Then we can add them 
to the virtual corrections to cancel the e poles, and ensure the virtual corrections are also 
integrable in four dimensions. This whole procedure can be illustrated by the formula 80| : 
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NLO 
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FIG. 5: Feynman diagrams for the massless quark emission contributions induced by SM QCD 
and NP interactions. 

where m is the number of final state particles at the LO, and da^ is a sum of the dipole 
terms. Besides, at hadron colliders, we have to include the well-known collinear subtraction 
counterterms in order to cancel the collinear divergences arising from the splitting processes 
of the initial state massless partons. Here we use the MS scheme and the corresponding 
NLO PDFs. 

For the process with two initial state hadrons, the dipole terms can be classified into four 
groups, the final-state emitter and final-state spectator type. 



1^ij,k{Pl, ■■■,Pm+l, 

1 



(Pi+Pj)^- 



mt 



i(..., tj, ..., k, ... 



Tu. ■ T,- 



■Vij^k\---,ij,---,k,...)m, (29) 



the final-state emitter and initial-state spectator type, 

T^tjiPl,--,Pni+i;Pa,--) = 



1 



.,a\---i "^3 }■■■■) O-i 



{Pi+Pjf -mf-Xij^a 

the initial-state emitter and final-state spectator type, 

Vf{pi,...,p^+i;Pa,...) = 

1 1 , ~ ~. ,T,-T 

-■{..., J,. ..■,m,... 



^^^^-^V^^.|...,^j,...;a,...)^,„ (30) 

ij 



^PaPi •^ij,a 



2 * J \---iJi ■••) "^5 -••I m,aii 



(31) 



10 



and the initial-state emitter and initial-state spectator type, 



P"*' {pi,...,Pm+i;Pa,Pb) 

1 1 



•^PaPi ■^i,ab 



~i{...;ai,b\ 



T,, ■ T, 



ai-irai,b 



\o-hb\ .„„• h\ ^ 



m,ai^ 



(32) 



where a, b and i,j, ... are the initial and final state partons, and T and V are the color charge 
operators and dipole functions acting on the LO amplitudes, respectively. The exphcit 
expressions for Xi^ab, Xij^a and V can be found in Ref. [80|]. The integrated dipole functions 
together with the collinear counterterms can be written in the following factorized form 



/ d^^"'\pa,Pb) m,ab{--; Pa, Pb\lm+a+b{£)\--; Pa, Pb)m,ab 

+ V / dx d<^^'^\xpa,Pb)m,a'b{--;XPa,Pb\P'^+bi^) +K.''^'ibi^)\- 
a' -^0 ^ 

+ {a ^ 6), 



XPa,Pb)m,a'b 



(33) 



where x is the momentum fraction of the splitting parton, d^^"^^ contains all the factors 
apart from the squared amplitudes, I, P, and K are insertion operators defined in 80|. 
For simplicity, in all the above formulas we do not show the jet functions that define the 
observables and are included in our numerical calculations. 

The operators P and K provide finite contributions to the NLO corrections, and only 
the operator I contains the IR divergences 



I 



IR 



as (47r)^ 
27rr(l-£) 



EEt.-t. 
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2 \ e 



1 



V{sjk, rrij, iTik] Sir) + --i2^j{mj, em) 



2 ( ^^ J V{sja, rrij, 0; Eir) + —^rj{mj,6iR) 

SjaJ J-j 



ry^Z 
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TleiR. 



-T„. ■ T, 



Sab 



2 \ e 



1 



1 7a 



with 



^IR T^a ^IR 



1 (Q]^' 



+ {a ^ b) 



(34) 



V{sjk,mj,mk;eiR) = — 

7i 



X ( 1 - 7TP, 



1 -2e 1 



2'^ 2 

Cf 



Pk 



2e 



-2 ' 
'IR 



T,{0,eiR) = ^, T^{m,^0,ejR) 

^IR ^IR 



(35) 



where Cp = 4/3, 7^ = 2, and ■jg = 11/2 — rif/S. And Sjk, Q'^j^, Vjk, and pn are kinematic 



11 



variables defined as follows 



Sjk = 2pjPk, Qjk = Sjk + mj +m^, Vjk = \ l- 



(PjPk) 



2' 



1 - Vjk + 2m2/(g2 _ ^2 ^ 

Pn = \ T-T. , o^2 ,rr^2 13:2 Z:2^ i^ = J,k). (36) 



^'-ml) 



1 + V,, + 2ml/ {Q% - m] - 



"^kJ 



When inserting Eq. flMj) into the LO amplitudes for the qq and qg subprocesses as shown 
in Eq. ( 133|) . we can see that the IR divergences, including the l/ejR terms, can be written 
as combinations of the LO color correlated squared amplitudes and all the IR divergences 
from the virtual corrections in Eq. (|26l) are canceled exactly, as we expected. 

IV. NUMERICAL RESULTS 

In the Lagrangian Cnp, there are totally nine free parameter Cli^, Cl^, C^l^ ^rr^ C'll; 
Clr, Cp^L, C^Yi ^^d A. If we include left-handed top quark t^ in the Cnp, it is suitable 
to work in SU{2)l doublet of the third-generation quarks (^l, ^l)"^- However, the Wilson 
coefficients Cll, Cll, C^ t an d C^-^ are highly constrained by the LEP data for the ratio of 
hb to hadron production 



56 = 0.121629 ± 0.00066, (37) 

which agrees well with the SM prediction. Thus, for simplicity we choose Cll = C*ll ~ 
^RL — ^KL — [ISl- Up to 0{a1/K^), contributions from color singlet operators due 
to mixing effects are much less than contributions from color octet operators, so we only 
consider color octet interactions. As a result, in the numerical calculations there are only 
three free NP parameters in the Lagrangian, i.e. Clr, Cp^j^ and A. 

Top quark mass is taken to be nit = 172.5 GeV. We choose CTEQ6L and CTEQ6M PDF 
sets ^j and the associated Ug functions for LO and NLO calculation, respectively. Both 
the renormalization and factorization scales are fixed to the top quark mass unless specified. 



We have used the modified MadDipole S^l package for the real corrections. 
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FIG. 6: Scale dependence of the total cross section at the Tevetron, the black and the red lines 
represent LO and NLO results, respectively. 

A. Scale dependence 

In Fig. [6] we show the scale dependence of the LO and NLO total cross section at the 
Tevatron for three cases: (1) the renormalization scale dependence fir = f^, f^f = ^^t, (2) the 
factorization scale dependence fir = ''^t. A*/ = Z^, and (3) total scale dependence fir = l^f = A*- 
It can be seen that the NLO corrections reduce the scale dependence significantly for all 
three cases, which makes the theoretical predictions more reliable. 

B. Tevatron constraints 



A-p-Q of top quark pair productions is defined as 
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FIG. 7: Values of C^^ and Cf^ allowed by the Tavetron data at la CL: (7(f=(7.50 it 0.48)pb 
and ^FB {mtt > 450 GeV)=0.475±0.114. The green star (5.31, -4.15) and red star (4.92, -3.88) 
represent the BFPs at LO and NLO level, respectively. 

where 



A 



NP 
FB 



/iSM 
^FB 



NP 



NP 



^NP ^NPn //^NP I ^ 

O^F -^B )/(^F +^: 



NP\ 



SM 



SM 



^F -^B )/(^F +^: 



.SM\ 



D ^NP//^bM I ^NF\ 



(38) 



are the asymmetries induced by NP and SM, and R is the fraction of NP contribution 
to the total cross section, a-p and o"b denote the total cross sections in the forward(F) 
and backward(B) rapidity regions, respectively. Up to order 0{a'j./A'^), total cross sections 
induced by NP can be written as 



a 



^E = [{OA28tr^l)Cl^ + {OA28l:'oZ)Cl^] (^)' pb, 



a 



_NP 

NLO 



[{OA42t^Zl)Cl,K + iOA35t'o'ol^)Cl^ (^)' Pb, 



(39) 



(40) 
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and the difference of the cross section in the forward and backward rapidity regions can be 
written as 

2 



CXt; 



NP _NP1 mtt>i50 GeV 
B JlO 



a 



[(o.ii8;°:°i)c|H,-(o.ii8+°:°i)c«i,](l^y pb, (41) 



[-1 



NP 



- (Tt 



Np-imtf>450 GeV 



B JnLO 



[(0.149+°;:^)C^^-(0.1031H^^)C^J (^^ pb, (42) 



where the errors are obtained by varying the scale between /i^ = /// = 'm(/2 and 
lir = fJ'f = 2mt. From the expressions Eqs. (139) - l42l) we can see that NLO corrections 
reduce the dependence of cxp^ — a^^ and a^^ on the renormahzation and factorization 
scales significantly. 

In Fig. [71 we show the allowed region in the (C^^r, Clr) plane that is consistent with the 
Tevatron data [9|: 



KX 



A 



EX 
FB 



(7.50±0.48)pb, 

0.475 ± 0.114, for ma > 450 GeV. 



(43) 



We use Monte Carlo programm MCFM [85] to get the cross section of the gluon fusion 
channel gg — )■ ti at the NLO QCD level. As for the process of qq — )■ ti, we have checked our 
value with the results given by MCFM at QCD NLO level, which are well consistent in the 
range of Monte Carlo integration error. Combining the contributions of these two channels 
we have the total cross section of ti production 



a. 



SM 



7 00+°-^^ 



pb. 



(44) 



where we have considered scale uncertainty in the calculations. For consistency we have 
used the SM predicted values of ApB imt > 450GeV) = 0.088 at NLO QCD level, although 
next-to-next-to-leading logarithmic (NNLL) SM QCD results are available |6|. In Fig. [TJ 
green and red regions correspond to NP LO and NLO results at la confidence level(CL), 
where we have considered theoretical and experimental uncertainty in the total cross section 
and only consider experimental uncertainty in the A-pB calculation. It can be seen that 
NLO corrections obviously change the allowed region of C^^^ and C^^. The green star 
(5.31,-4.15) and red star (4.92,-3.88) represent the best-fit point(BFP) at LO and NLO 
level, respectively, from which we can see higher order corrections reduce the BFP by about 
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7%. The tt total cross sections induced by NP at the NLO BFP (4.92, -3.88) are 



NP 



cT^fLo = 0.445 pb, 



<NLO = 0.497 pb, 



(45) 



where the K factor is about 1.12. Ap^ containing NP contributions at the NLO BFP are 
shown together in Table HI The NLO QCD corrections to A-p^ can reach about 10%, and 
the theoretical predictions containing NP NLO effects are consistent with experimental 
results at 2a CL. 

In Fig. [HI we show differential cross section da/drriu when we consider NP effects at 





SM NLO QCD + NP LO 


SM NLO QCD + NP NLO 


app 


0.175 


0.189 (~ 0.7 cr) 


Ati 


0.252 


0.275 (~ 1.6 cr) 


ApB i'^tt < 450 GeV) 


0.132 


0.136 (~ 1.6 cr) 


^FB i^tt > 450 GeV) 


0.452 


0.475 (~ a) 


^«B (|Ay| < 1) 


0.170 


0.161 (~ 0.9 cr) 


^«B {\Ay\ > 1) 


0.719 


0.681 (~ 0.3 cr) 



TABLE L Afb with C|p(lTeV/A)2 = 4.92 and Cl^{lTeV/Af = -3.88 at the Tevtron, where 
^pg and Ap-Q are the ApB in the lab frame and the tt rest frame, respectively, and Ay = yt — Vt 
is the difference of rapidities of the top and antitop quarks. Here we list the CL when containing 
NP effects at NLO level. 

the NLO BFP, from which we can see that higher order corrections do not change the 
distribution very much. 



C. LHC predictions 



The process of top quark pair production has been measured at the LHC, and the cross 



section [87, 



IS 



af^^^{V^ = 7 TeV) = 180 ± 18pb, 
a'i^^iV^ = 7 TeV) = 158 ± 19pb, 



(46) 
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FIG. 8: Differential cross sections da /drriii as a function of rriii at the NLO BFP (4.92, —3.88). 
Here " Exjaerimental data" is da/dnif-i measured with 2.7 fb^"^ of integrated luminosity at the 
Tevatron (86']. "SM NLO QCD" represents the results in the SM QCD at NLO level. "NP LO + 
SM NLO QCD" and "NP NLO + SM NLO QCD" stand for the predictions including NP effects 
up to 0{as/ h?) and ©(a^/A^), respectively. 



which is consistent with the SM predictions. The NP contributions at the NLO BFP 
(4.92, —3.88) is about 3 pb, which is much smaller than the experimental uncertainty. Thus, 
it is difficult to measure the NP effects only through the cross section measurements. In 
Fig. [HI we show differential cross section da/dmu at the LHC when we consider NP effects 
at the NLO BFP (4.92, —3.88), from which we can see that NP contributions almost do not 
change the distribution. 

Since the LHC is a proton-proton collider, which is forward-backward symmetric, the Ap^ 
defined at Tevatron can not be directly applied to the proton-proton collider experiments 



at the LHC. The Ac used by CMS |89|] can be written as 

^ ^ (^i\vt\ - \m\ > 0) - (^{\vt\ - \vi\ < 0) 

^ a{\7]t\ - \7]t\ > 0) + a{\7]t\ - \7]i\ < 0) ' 



(47) 
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FIG. 9: Differential cross sections da/dm^f as a function of m^f at the LHC with ^/S = 7 TeV. 

where rjt and rjt are pseudo rapidities of top and antitop quark, respectively. Its value is 
measured to be 



Ac = -0.016 ± 0.030(stat.)+°;°i°(syst.) 



(48) 



which is consistent with the SM predictions: Aq = 0.013(11) [l7|, l89|, l90|] . The Aq induced 
by NP interactions at the NP NLO BFP (4.92, -3.88) is 0.063, which is about 5 times larger 
than SM predictions, and also consistent with the CMS data at about 2a CL. 

In Fig. [ini we show the results of a combined fit to the ti data in the presence of NP at 
different CLs. The blue contours from dark to light indicate the experimentally preferred 
region of 68%, 95% and 99% CL in the (C^^, C|^) plane. The black dot represent the SM 
point (0,0), and the black star represent the NP NLO BFP (4.92, -3.88). The black doted 
and dashed lines respectively represent the value of ti cross section and the Aq at the LHC 
with y/S = 7 TeV. From the location of the blue area, one finds that Ac predicted by NP 
is obviously much larger than SM predictions, and LHC has potential to find this difference 
when the measurement precision increases. 
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FIG. 10: Results of a combined fit to a^i and the value of ^fb allowing for NP at different 
CLs. The blue contours from dark to light indicate the experimentally favored region of 68%, 
95% and 99% CL in the (C'^iJ'^f/?) pl^-iie. The black doted and dashed lines respectively 
represent the value of tt cross section and the Aq at the LHC with \fS = 7 TeV. The black 
dot and black star represent the SM point and the NP NLO BFP. 

V. CONCLUSIONS 

In conclusion, we have investigated ApB and total cross section of top quark pair pro- 
duction induced by dimension- six four quark operators at the Tevatron up to 0{a1/A'^). 
Our results show that, NLO QCD corrections can change ApB and the total cross section 
by about 10%. Moreover, NLO QCD corrections reduce the dependence of ApB and the 
total cross sections on the renormalization and factorization scales significantly, which lead 
to increased confidence on the theoretical predictions. We also evaluate total cross section 
and Ac induced by these operators at the LHC up to (9(a^/A^), for the parameter space 
allowed by the Tevatron data. We find that the value of Aq induced by these operators is 



19 



much larger than SM predictions, and LHC has potential to discover these NP effects when 
the measurement precision increases. 
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Appendix A: FINITE TERMS IN VIRTUAL CORRECTIONS 

In this appendix, we collect explicit expressions of finite terms in virtual corrections. 
Analytical continuation for the Mandelstam variables are defined as 

s — )> s + ie, 
u — )■ u + ie, 
t -> t + ie. (Al) 

For simplicity, we introduce the following abbreviations. 
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First we list the amplitude for helicity (H 1-+), 

<NP,fin/ , , , ^ iasin 

+Mos{+-++)Ci+Moo{+-++)Cs\, (A7) 



zA^::;^^ °(+ - ++) = -^Cp [Mssi+ - ++)Ci + MU+ - ++)Cs 



where 



+Crr [h9{s)Mi + hs{s)M2 + 8h{s) (A^/™) + A^2^"^^)] , (A8) 

-Mso(+-++) = -Clj^[f2iu)Mi + Mu,t)M2- fiiu) {M,^"'^ + M2^"'^)] 

+Crr [Mt,u)Mi + f2it)M2 - flit) (A^/™) + A^2(™))] , (A9) 

-Mos(+-++) = -lcl^[f2iu)Mi + Uu,t)M2-fiiu){Mi^"^^+M2^"'^)] 

+lcL [Mt)Mi + f2{t)M2 - flit) (A^i^"^) + M2^"'^)] , (AlO) 

-Moo(+ - ++) = C^L [gsis, u)Mi + gi (s, ^(t), /3H) M2 + ^5(5, u) (A^i^™) + ^^2^"^^)] 
+CIk [91 is, fsit), f,iu)) Ml + g2is, t)M2 + ^4(5, t) (A^/™) + ^^2^'"^)] 
+Cl^his)Mi + C^l/^6(s)-M2. (All) 

Similarly, amplitude for helicity ( — h ++) can be written as 
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Appendix B: evolution equation of Wilson coefficients 

In this appendix we present the evaluation equations of the Wilson coefficients C^g and 
C|^g, expanded to 0{as)- 

Cldf^r) = C^lK) - 0.0397887C^a,ln^C«L(m,), (Bl) 
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